Superfluid Hydrodynamics of an Electron Gas in a Superstrong Magnetic Field 
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We derive the equations of hydrodynamics of a fully polarized electron gas placed in a strong 
magnetic field. These equations reveal the existence of solitons - immobile or propagating domain 
wall-like defects whose plane is perpendicular to the field direction. The solitons are used to construct 
weakly excited stationary states, and novel nonuniform persistent current states of the system. 
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Understanding the properties of matter subject to an 
extremely strong magnetic field Q has fundamental im- 
portance. At fields of order IQ^G the interaction of 
atomic electrons with the external magnetic field begins 
to dominate over their Coulomb interaction which results 
in a novel kind of atomic physics In semiconduc- 
tors and semimetals due to the smallness of the effective 
mass of the charge carriers and dielectric screening of the 
Coulomb interaction the corresponding field is sub- 
stantially smaller, and accessible in the laboratory [Q. 

Since a uniform magnetic field localizes the motion 
of a charged particle in the plane perpendicular to the 
field direction, and does not affect the motion along 
the field ||^, strong fields can be used to study essen- 
tially one-dimensional physics using three-dimensional 
samples. The one- dimensionality is most pronounced if 
the field is strong enough to pull all the carriers into the 
lowest Landau level. In this limit in an interacting sys- 
tem a variety of phases of one-dimensional nature such as 
charge-density waves , Wigner crystals ^ , and ex- 

citonic insulators 1^ have been predicted, and a uniform- 
density state (adiabatically evolving from the noninter- 
acting limit) occupies only a portion of the phase diagram 

. Recent work has demonstrated that this phase has 
Luttinger- liquid-like properties . 

In this Letter, by studying elementary excitations, we 
show that the uniform-density phase behaves as a uni- 
directional superfluid. The governing equations of hy- 
drodynamics have a rich mathematical structure which 
owes its existence to the one-dimensional physics of the 
problem. A detailed treatment will be presented else- 
where [|lO| . Below we restrict ourselves only to the salient 
physics effects. 

We start with reminding the reader of the method of 
constructing the ground state of a noninteracting elec- 
tron gas of density n placed in a fully polarizing magnetic 
field H pointing along the x direction of a coordinate sys- 
tem. The motion of electrons is confined within infinitely 
long Landau tubes parallel to the field direction: the mo- 
tion perpendicular to the tube axis is localized within the 
magnetic length I = {hc/eHy^^ while the motion along 
the tube is free. The degeneracy of the lowest Landau 
level, i.e. the number of Landau tubes in the system of 
cross-sectional area S perpendicular to the field direction 
is given by § eHS/2TTfic = S/2t:P which imphes that 



the cross-sectional area of the individual tube is 27rP. 
Therefore the magnetic flux confined within the tube is 
the fiux quantum he/ e, and thus each Landau tube can 
be thought of as a magnetic field line. Hence the quantum 
state of every electron can be characterized by specify- 
ing the Landau tube it belongs to, and the wavevector 
kx describing the free motion along the field. Because 
of the Pauli principle, for a given Landau tube all the 
states with \kx\ < kp are occupied while all the states 
having | fc^; | > kp are empty. The Fermi wave vector k p 
can be found as kp = np where p is the linear electron 
density in the tube, which can be computed by divid- 
ing the total number of electrons N by the number of 
Landau tubes S/2ttP, and by the tube length L. Thus 
p — 2'iTpn, kp = 27r^Z^n and the Fermi energy can be 
found as ep — h^kj^/2m — {2n^hPn)'^ /2m where m is 
the electron mass. The expression for the ground-state 
energy 
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(where V — SL is the system volume) can be rewritten 
as (L7r^?i^/6m)(27rPn)^ times S/2nl'^: the ground-state 
energy of a one-dimensional spinless Fermi gas of linear 
density p = 2'!rl'^n, multiplied by the total number of Lan- 
dau tubes. The condition of being at the lowest Landau 
level is satisfied if the Fermi energy lies below the bot- 
tom of the second Landau level fieH/mc [|| . This leads 
to the condition nl^ < l/7r^\/2 which gives the range of 
applicability of Eq.(|l]). In what follows we will adopt 
the stronger condition nl'^ <C 1 which also allows us to 
neglect transitions to the higher Landau levels. 

The presence of the Coulomb interaction introduces a 
new length scale into the problem, the "exciton" Bohr 
radius a — eh^ /me^, where e is the dielectric constant of 
the medium. In astrophysical applications when e = 1, 
and m — m^, the electron mass in vacuum, a reduces 
to Bohr's radius. In the limit nl^ ^ 1 the average dis- 
tance between the electrons belonging to the same Lan- 
dau tube, l/2iTpn, is much bigger than the average dis- 
tance between the nearest tubes I. Thus the electrons 
belonging to the same Landau tube interact much more 
weakly among themselves than the electrons belonging 
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to different tubes. Therefore the physical picture of in- 
dependent Landau tubes survives if the Fermi energy is 
much bigger than the Coulomb interaction of electrons 
confined to different tubes. The latter (in an electrically 
neutral system) is estimated as the plasmon en- 

ergy hujp where uj^ = ATrne^/meQ is the plasma frequency. 
Combining the requirement S> hiOp with the condition 
that all electrons are at the bottom of the lowest Landau 
level, we arrive at the concentration interval 
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which is the range of applicability of our theory. The 
interval (^ exists if the exciton Bohr radius a is much 
bigger than the magnetic length I; the magnetic fields 
satisfying this condition are refered to as superstrong M . 
Both length scales coincide at the field Hq = m^e^c/e^fi . 
In astrophysical applications (m ~ mg, e = 1) = 2.35- 
lO^G, while in condensed matter applications the corre- 
sponding field is {eme/rnf ^ 10^ — 10^ times smaller, 
and for a/Z ~ 5 — 10 the density range (|^) exists in doped 
semiconductors and semimetals. 

We postulate that the long- wavelength physics of a po- 
larized electron gas placed in a slowly varying external 
potential V{y) is described by the energy functional of 
the complex order parameter ^'(r, t): 
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The physical meaning of ^' is that = n(r, t) gives the 
electron density which is a function of position r and time 
t. If ^! is restricted to be real, then (^) can be recognized 
as a starting point of a generalized Thomas-Fermi density 
functional theory ||l^: for the translationally-invariant 
case the density is uniform, and then (H) reduces to the 
ground-state energy (|l|). The second term of (^ is the 
energy in the external field. The derivative term, known 
as "Weizsacker inhomogeneity correction" [Q, which is 
a lowest-order term of a gradient expansion, takes into 
account the fact that deviations from uniformity in the 
field direction are costly. Density functionals similar to 
(^ are used to describe static properties of multielcctron 
atoms placed in very strong magnetic fields ||]. 

We further postulate that the order parameter ^ 
evolves according to the equation of motion 
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In the representation 5" 
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Eq.(W) reduces to the system of two coupled equations: 



dtn + dx{nv) = 
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In the long-wavelength limit when the last term in the 
right-hand side of (||) can be dropped, the system of equa- 
tions (H) and (^) can be identified with those of hydrody- 
namics of an ideal fluid whose pressure and density 
are related by the equation of state p = {2-K'^hl'^)'^n^ /2>m. 
The variable v = {h/m)dxO has a meaning of the super- 
fluid velocity which only has a component along the field 
direction. We note that the equations of hydrodynamics 
of an electron gas in a superstrong magnetic field could 
be written down from the outset and the density depen- 
dence of pressure can be computed from the ground-state 
energy (|l]) a,s p — —dE/dV. This provides an important 
test of the underlying equations (||) and (^); specifically, 
selection of the gradient term of (||) in the Weizsacker 
form plays a central role in getting the correct hydrody- 
namical limit. We stress that hydrodynamical descrip- 
tion requires two fields, density and velocity; that is why 
the order parameter must be complex. 

Spectrum of small density oscillations: - Let us look 
for solutions to (^ (with = 0) of the form 'I'(x,t) = 
ip^x, t) exp(—ifit/h). The function ')p{x, t) then obeys the 
equation 
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+ (2^2/2f(|V'|^-Vo)]V' (7) 



where the chemical potential /i = T:^h'^{2'Kp)'^^Q/2m 
is selected so that the ground state corresponds to the 
uniform electron density no = V'o- Consider a small 
oscillation of 'ijj around its constant value ipo- ip = 
'ipo + Aeiq)[i{qx — LL!t)] + B* e:cp[—i{qx — ut)], where A and 
B* are small complex amplitudes, and q is the wavevector 
along the field direction. Substituting the expression for 
ip into (0), linearizing, and solving the resulting system 
of two equations, we find for the dispersion law 
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The spectrum takes the Bogoliubov form [Q: for large 
momenta the dispersion law coincides with the free- 
particle energy h^q^ /2m, while at low momenta we get 
instead the phonon dispersion lo = eg, where c = 
2TT'^PfinQ/m is the sound velocity. We note that this 
velocity coincides with the hydrodynamical expression 
c = (c?p/m(in)i/2 for a gas with equation of state p — 
{27r'^hP)'^n^ /3m. The Landau criterion applied to (||) 
then implies that the critical velocity of superfluid 
motion along the direction of magnetic field equals the 
sound velocity. 

Solitons: - Without the external potential {V — 0) 
Eq.(^ has soliton solutions These describe moving 
(gray soliton) or standing (dark soliton) density depres- 
sions accompanied by a corresponding profile of phase 
p6[ . Because of the system uniformity in the plane per- 
pendicular to the magnetic field, the solitons are domain 
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wall-like defects whose plane is perpendicular to the di- 
rection of magnetic field; the solitons can only propagate 
along (or opposite) to the field direction (Figure 1). 

n 




FIG. 1. Schematic illustration of the density profile due to 
the soliton: electron density depends only on the coordinate 
X in the direction of the magnetic field H\ transverse to the 
field direction is labeled (y, z). The phase profile (not shown) 
is a smeared step function (|lO[); the step position coincides 
with that of the density dip. 

In dimcnsionless variables v — u/hq, y = 27r^/^noa;, 
r = (27r^Pno)^?ii/m, C, — y — /3t, the density {v) and 
phase {9) profiles for the soliton are given by 
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where the dimcnsionless velocity (3 is measured in units 
of the sound velocity. These formulas imply that as 
(3 varies between zero and unity, the soliton becomes 
more shallow, and the phase shift across the soliton 
AO = cos-i[(3/32 - 1)/(1 + 3/32)1/2] changes from -tt 
(/3 = 0) to zero (/3 = 1); the soliton solution ceases to ex- 
ist for velocities exceeding the sound velocity. Microscop- 
ically the domain wall is composed of aligned point-like 
solitons belonging to different Landau tubes, and there is 
a positive energy per unit area associated with the planar 
defect. It can be calculated as the energy of the soliton 
belonging to the Landau tube (computed in ||T^), divided 
by the cross-sectional area of the Landau tube 27rZ2 : 



a(/5) = ^^^^i^(l-/32)ln 
m 
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The positivity of the surface energy ( [TTI ) guarantees the 
stability of the domain wall (|9|), ( |l0|) against transverse 
fluctuations. 

Stationary states: - Solitons form a natural language 
to describe weakly excited states of the system. For ex- 
ample, let us look at electron gas placed inside a layer 
of macroscopic thickness L; the plane of the layer is per- 
pendicular to the direction of the magnetic field. If the 
electrons cannot leave the layer, then in the ground state 
the density profile will look as shown schematically in 



Figure 2 (left): the density stays constant almost every- 
where except for the narrow vicinity of the walls where 
it sharply plunges to zero. 
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FIG. 2. Sketch of the first (left) and the second (right) 
excited states of electron gas in a box. 

It is convenient to think of these two dips as having 
one dark soliton present in the system. It is also con- 
venient to take the convention that what is shown in 
Figure 2 (left) is a first excited state with respect to the 
strictly uniform state which we will call the ground state. 
Thus the first excited state has just one soliton present 
in the system. Then the second excited state will have 
two solitons [Figure 2 (right)]. The argument can be ex- 
tended to any number of solitons. For example, if there 
are Q dark solitons present, then as long as soliton tails 
do not overlap, Q <^ PnoL, the system energy is given by 
E{Q) — Eq + S<t{Q)Q where Eq is the ground-state en- 
ergy (|^) . A similar conclusion is applicable to the case of 
a ring-shaped sample with the only constraint that now 
Q must be even because the phase shift due to the dark 
soliton is ±7r. Having an even number of dark solitons 
will guarantee that the overall change of phase of the or- 
der parameter upon going around the ring is multiple of 
27r, and thus the order parameter is unique. 

Because translational symmetry is not broken for the 
ring geometry, dark soliton states do not exhaust all pos- 
sible excited states of the system: gray solitons can also 
be involved. However it will be easier to anticipate their 
role after discussing the persistent current states. 

Persistent currents: - Persistent currents are only pos- 
sible in the ring geometry. Experimentally the magnetic 
field directed along the circumference of a ring-shaped 
sample can be generated by an electric current flowing 
along the symmetry axis perpendicular to the plane of 
the ring. In the ring geometry the phase of the order 
parameter can only change by a multiple of 27r. The sim- 
plest way to accommodate a nonzero phase change is to 
have a constant phase gradient along the ring [0 . This 
will lead to a spatially uniform persistent current state in 
which all the electrons are moving around the ring with 
the same velocity. This state can be also viewed as a 
stationary state. 

There is also a novel nonuniform way of accommodat- 
ing the phase difference. The corresponding persistent 
current states can be constructed from gray solitons, and 
they have a series of remarkable properties. For exam- 
ple, imagine we want to insert into the ring Q solitons 
to accommodate an overall phase shift of 27r. Then the 
phase shift due to the individual soliton must be 2tt/Q. 
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This will in turn determine the soliton velocity (3 which 
will then fix the surface energy ( p] ) . An example of four 
gray solitons moving along the ring is sketched in Figure 
3. 




FIG. 3. Schematic illustration of four solitons traveling on 
a ring. The upper curve is the density profile, and the lower 
curve indicates the zero-density level. The arrow shows the 
direction of motion of the soliton train. 

The energy of the train consisting of Q solitons run- 
ning around the ring of circumference L is given by 
E{Q) — Eq + S(j[l3{Q)\Q as long as soliton tails do not 
overlap, Q <C (1 — f3'^y^''^l'^noL. In the reference frame 
moving with the soliton train the solitons will be stand- 
ing but the background electron gas will be moving with 
a constant velocity. This gives an example of a stationary 
state involving gray solitons. 

For the uniform moving state the total angular mo- 
mentum of the system is quantized in units of Nh where 
N is the total particle number [0. We speculate that 
for the soliton train moving state the total angular mo- 
mentum might be a fraction of Nh. Our contention is 
based on the observation that while the soliton is mov- 
ing, the electron motion is localized only in the vicinity 
of the soliton core - merely a tiny fraction of particles are 
involved in the motion. 

In applying these results to the real system one has to 
keep in mind that since a soliton is a density depletion, 
for solitons which do not travel fast enough, the den- 
sity at the soliton core can drop below the left bound of 
. Therefore our results are strictly applicable only to 
solitons whose velocity is larger than some critical veloc- 
ity. Future work is necessary to understand how interac- 
tions between the electrons belonging to different Landau 
tubes affect the slow solitons. One possibility is that they 
convert the depleted soliton core into a two-dimensional 
Wigner crystal. Another possibility is a modulational 
instability of the soliton in which Coulomb repulsion of 
the electrons decreases at the expense of increase of the 
surface energy (|ll|). 

The unique mathematical structure of the underlying 
equations (^-(^ which so far did not enter the discussion 
deserves a special commentary. First we note that in the 
hydrodynamical limit when the last term of can be 
neglected (and = 0) the problem of arbitrary motion 
of a one-dimensional gas with the equation of state p ^ 



can be solved in a closed form jisj. Second, (even 
for V ^ x'^) equations (U)-(||) have self-similar solutions 
]l9t which only exist for the equation of state p ^ n^. 
Physically some of them describe free expansion of the 
electron gas initially confined along the field direction 
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